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^^ . Abstract. We study the generation of helical magnetic fields during single field 

r) I inflation induced by an axial coupling of the electromagnetic field to the inflaton. 

^ ' During slow roll inflation, we find that such a coupling always leads to a blue spectrum 

Q^. with B^{k) oc fc, as long as the theory is treated perturbatively. The magnetic energy 
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density at the end of inflation is found to be typically too small to backreact on 



^ ■ the background dynamics of the inflaton. We also show that a short deviation from 

C/5 , slow roll docs not result in strong modifications to the shape of the spectrum. We 

^ ' calculate the evolution of the correlation length and the field amplitude during the 

inverse cascade and viscous damping of the helical magnetic field in the radiation 

^^ \ era after inflation. We conclude that except for low scale inflation with very strong 

p^ ' coupling, the magnetic fields generated by such an axial coupling in single field slow 

^sj , roll inflation with perturbative coupling to the inflaton arc too weak to provide the 

^^ ' seeds for the observed flelds in galaxies and clusters. 
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1. Introduction 

Cosmic magnetic fields liave been observed on all scales ranging from stars to near and 
;-i ■ far away galaxies and galaxy clusters [H El [3l Hj. The strength of the magnetic fields 

observed in galaxies and clusters is typically of the order of /iGauss. Recently, using 
the absence of extended GeV emission around TeV blazar gamma-rays, a lower limit of 
3 X 10'^® Gauss on the strength of intergalactic magnetic fields was derived [5], El [71 E] ■ 
These observations prompt the question of the origin of cosmic magnetic fields. 
Have they been generated during structure formation or have primordial magnetic 
fields been amplified? So far, this question has no clear, satisfactory answer. Several 
studies [HI [ini [m [121 [131 [HI [iSl [IE] of second order perturbation theory have shown 
that up to recombination only very weak magnetic fields of the order of lO"^*^ Gauss 
can be generated by structure formation. However, it turns out that the magnetic fields 
from second order perturbation theory are not strong enough to exceed the lower limit 
derived in Refs. [3 [6l [71 [8] . Whether such fields could have been generated later by the 
process of galaxy formation and then ejected into intergalactic space remains unclear. 



In this work, we pursue the idea that instead magnetic fields are of primordial 
origin and might have been generated in the early universe. Primordial magnetic fields 
are interesting as they induce all three kinds of gravitational perturbations i.e. scalar, 
vector and tensor; all of which contribute to the Cosmic Microwave Background (CMB) 
temperature and polarisation anisotropies. These primordial fields also lead to non- 
Gaussian signals in the CMB even at the lowest order contrary to the higher order 
effect due to infiationary scalar perturbations. Analyses of such effects using recent 
cosmological data provide an upper limit of a few nCauss on the primordial magnetic 
fields [171 Ell UHl ISO]- It has been argued that the electroweak or the QCD phase 
transitions, if they are first order, lead to the generation of magnetic fields [211 I22l[23ll2^ . 
Also non-perturbative processes during preheating can generate significant magnetic 
fields with, in some cases, a helical component, see for instance [25]. However, causality 
strongly constrains such fields. Their power spectrum is very blue with B^{k) oc k"^, and 
therefore, their amplitude on large scales is far too small [261 127] . 

On the other hand, if the magnetic fields are produced during infiation, their power 
spectrum is a priori not constrained by causality but only by the specific model. Since 
the standard electromagnetic (EM) action is conformally invariant the fiuctuations in 
the EM field are not amplified in the conformally flat expanding background of infiation. 
In order to generate magnetic fields, one needs to break conformal invariance of the EM 
field, e.g. by coupling the EM field to a scalar or a pseudo-scalar field or to a curvature 
invariant (for an overview see, for instance, Ref. [28]). Typically, a term of the form fF"^ 
is considered, where / is a function of time or of the infiaton and F is the EM field tensor. 
Depending on the form of the coupling /, this gives rise to different magnetic field power 
spectra, and even scale-invariant spectra are possible (in this context, see Refs. [29ll30] ). 
In this way, the magnetic fields can have sufficient amplitude on large scales to provide 
seeds for the observed fields in galaxies and clusters. However, the backreaction on the 
infiation dynamics, the production of gravitational waves and nucleosynthesis bounds 
on the amplitude of gravitational waves strongly constrain the magnetic energy density 

[niEZ]. 

Quite a different situation is encountered if a coupling to the parity-violating term 
FF, i.e. a term fFF, is added to the standard EM action F^, where F is the dual of F. 
As a consequence, magnetic helicity is generated, which is absent in the case discussed 
above. This has two interesting consequences: firstly, contrary to non-helical fields, 
helical fields evolve in the cosmic magnetohydrodynamic (MHD) plasma via inverse 
cascade [311 |32l I33] . This transfers power from small to large scales so that even blue 
spectra can lead to significant power on large scales. In Ref. [34], it was shown that the 
inverse cascade is not quite sufficient for helical fields generated at the electroweak phase 
transition [35], but it might work for magnetic fields from infiation. Secondly, helical 
magnetic fields leave a very distinct signature as they violate parity symmetry. This 
leads to observable effects, e.g. correlations between the anisotropies in the temperature 
and B-polarisation or in the E- and the B-polarisations in the CMB [36]. Furthermore, 
they induce helical gravitational waves [36] which might be observable 



Some consequences of primordial helical magnetic fields and their generation from 
primordial helicity have been studied in the past [381 EH]- The interactions of helical 
magnetic fields and axions have also been investigated jlQlHT]. Recently, the generation 
of helical magnetic fields during inflation in specific models has been studied, e.g. in 
Ref . |12] , helical magnetic fields from N-fiation were investigated. In this case, the large 
number of pseudo-scalar fields driving inflation effectively leads to a large coupling, 
/ oc viV, to the FF term. In Ref. |13], some toy models for the coupling were analysed 
where / was taken to be a power law function of kt [k being the wave number and t the 
conformal time). 

In this paper, we study magnetic fields generated by an axial coupling of the form 
f{(f))FF during inflation, where (p is the infiaton. We consider two different forms of the 
coupling function and show that, contrary to a non-helical coupling of the form /(0)F^, 
a helical coupling always leads to a spectral index n = 1 for B'^{k) oc /c", as long as slow 
roll inflation is considered. We derive the condition for the theory to be perturbative, 
i.e. the free part of the action dominates over the interaction. Of course this is not 
necessarily true, but at least naively if not, we can no longer trust our calculation of 
particle creation out of the quantum vacuum which is based on perturbative quantum 
field theory. We estimate the magnetic energy density as a function of scale and show 
that backreaction is typically small. These conclusions are valid for any reasonable 
coupling function /(0). We confirm our analytical results numerically. Furthermore, 
we study the effects of a short deviation from slow roll on the magnetic field spectrum 
and show that such deviations, if kept within the bounds permitted by the CMB data, 
do not strongly modify its shape. Even though the inverse cascade in the radiation 
dominated era after inflation does move power to larger scales, the final strength of the 
magnetic field on cosmologically interesting scales is still insufficient to provide seeds for 
the observed magnetic fields in galaxies and clusters, except if the inflation scale is low, 
T^, < lO^GeV and the axial coupling is very strong. Even if we assume very efficient 
dynamo amplification. 

The remainder of this paper is organised as follows. In the next section, we 
introduce the axial coupling of the EM field to the infiaton, derive the field equations, 
discuss the background evolution and the slow roll approximation and derive the linear 
perturbation equations of the infiaton and the EM vector potential. In Sec. [3l we discuss 
the evolution equation of EM quantum fluctuations during inflation and compare the 
helical to the normal case. In Sec. HI we derive a condition on the chiral coupling by 
requiring the theory to be perturbative, investigate different coupling functions in the 
slow roll approximation, and solve the evolution of the EM fluctuations analytically. 
In Sec. [5l we discuss the consequences of a brief violation of the slow roll approximation 
on the magnetic field spectrum. In Sec. El we study the evolution of the magnetic 
field during the radiation era after inflation and determine the final spectrum after the 
inverse cascade. Finally, in Sec. [71 we conclude with a few comments on our results. 
Three appendices contain some details on the quantisation of the vector potential, 
the condition on the coupling function from the perturbativity of the theory and the 



asymptotic behaviour of the Coulomb wave functions, respectively. 

Notation and units: We work in a metric with signature ( h + +). For tensor 

components, Greek indices take values ... 3, while Latin indices run from 1 to 3. 
The components of spatial 3-vectors with respect to a comoving basis are denoted in 
bold face. We employ Heaviside-Lorentz units such that c = h = ks = e^ = fio = 1. 
The reduced Planck mass is defined as nip = (SvrG)"^/^. We normalise the cosmic scale 
factor to unity today so that the comoving scales become physical scales today. 

2. Axial coupling of electromagnetism to the infiaton 

2.1. Action and field equations 

We consider a scalar field, 0, which takes the role of the infiaton and the EM field, 
F^y = df^Ai, — dyA^, characterised by its four- vector potential A^. The EM field is 
conformally coupled to the metric and therefore, no fluctuations are generated unless 
there is either an explicit coupling to the infiaton or conformal symmetry is broken 
directly, e.g. by coupling F to a curvature term. Here, we investigate the first possibility 
and study a helical couphng given by the action 

S[ct>.A^] = fd'x,/^{CM + ^UA,) + Cii<p,A^)}. (1) 

The Lagrangian densities of the free fields are 

C^{<P) ^ ^g-0(d^^)(dp^) + V{<P) (2) 

CeUA^) = -^Fa/3F°^ (3) 

and the axial interaction is given by 

Cj{<f>,A,) ^l/(0)F„^F"^ (4) 

It describes a coupling of the scalar field to the parity-violating term, FF, where F is 
the dual of the EM field tensor and is defined as 

F'^" = Iv'^'^^F^^- (5) 

Here 77^'^"'^ is the totally anti-symmetric tensor in four dimensions with /^o^^s = (■_(^^-i/2_ 
For an observer with 4-velocity u'^ the electric and magnetic fields are F^ = F^oiU°' and 
B^ = Ff^aU°', respectively, and we have F^iiF"^ = —AEaB'^. 

If the scale of inflation is above the electroweak scale then, in principle, one should 
specify whether (j) couples to f/(l), SU{2) or both of the electroweak SU{2)l x f/(l)y. If 
the coupling is to both, which seems simplest, the same process that leads to magnetic 
field helicity also induces a non-zero baryon number which is related to the electroweak 
Chern-Simons number [35]. The electromagnetic Chern-Simons number is equivalent to 
the hehcity. In this work we do not discuss this additional potentially interesting aspect, 
but concentrate on the magnetic fields remaining after the end of the electroweak phase 



transition. We assume the conversion from W^ and Y into photons to be efficient and 
not to affect the resulting magnetic field distribution significantly so that we may simply 
consider the coupling of to the photon field. 

The axial coupling is characterised by the scalar function /(0). We will see later 
how this function affects the evolution of the vector potential. Notice that if /(</>) was a 
constant, the EM part of the action would still be conformally invariant and therefore, 
no EM fluctuations could be amplified during inflation. Note also that the term £/ 
either breaks parity explicitly if is a normal scalar field or, if is a pseudo-scalar, 
parity is broken spontaneously by the presence of a background field 0^0. For the 
discussion in this work, this distinction is not relevant. However, in certain models, it 
might be relevant for the amount of parity violation generated during reheating. 

Varying the action with respect to leads to a sourced equation of motion for the 
scalar field 

V"9„0-\/'(0) = ^/'(0)F«;3F"^ (6) 

The primes in V'{(j)) and /'(0) denote derivatives with respect to 0. The field equations 
for the EM field follow from varying the action with respect to A^: 

V„F^" = /'(0)a0)F'^". (7) 

Comparison with the usual inhomogeneous Maxwell equation leads us to interpret 
the source term on the right hand side as an effective axial currently. To obtain the 
above form of the inhomogeneous Maxwell equation, we used the homogeneous Maxwell 
equation 

V^F'^" = ^ V[aFh = (8) 

which is equivalent to the Bianchi identity, dF = 0. 

2.2. Background evolution 

To describe the universe during inflation, we work in a flat Friedmann-Lemaitre (FL) 
background metric characterised by the line element 

rfs^ = a^{t){-df + 6ijdx'dx^) (9) 

where a is the scale factor and t is conformal time which is related to cosmic time r 
hj adt = dr. Derivatives with respect to conformal time are denoted by a dot and the 
conformal Hubble parameter is "H = a/ a = aH where H = {da/dT)/a is the physical 
Hubble parameter. 

We assume the scalar field to dominate the energy budget of the universe and to 
drive inflation. We shall check later under what conditions the EM energy density is 
negligible and this approach is justified. We decompose the scalar field into a background 
value and a small perturbation: 0(x'^) = ip(t) + 5(j){x'^). At background level, Eq. (P) 
reduces to the homogeneous evolution equation for the scalar field 

(p + 2n^ + a^V'{^) = 0. (10) 

I An axial anomaly which also induces a source term of this form has been discussed in Ref . |44j . 
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The evolution of the scale factor is determined by the background Friedmann constraint 
equation 

3mln^ = ^^^ + a^V{^). (11) 

During slow roll inflation, the potential of the inflaton fleld is dominating the energy 
density of the universe and the flrst term on the right hand side of Eq. (TTT]) can be 
neglected. This is quantifled by means of the slow roll parameters (see e.g. [15]) 



mi fV . 



' ml V" 



3 V « ^^ (^2' 

To flrst order in the slow roll parameters, the evolution equation of (f and the Friedmann 
constraint can be reduced to [45] 

V9 ^ -^ = ±V2'e£^ ^ ±V2'emp'H. (13) 

The sign of if depends on the details of the inflation model, but here we do not assume 
a speciflc form of the potential. Using also that "H ~ |t|~^ and that e is roughly constant 
in the slow roll regime, we can integrate this result to flnd 

if ^ v9in=F V2impln(t/tin) (14) 

where ipi^ = '^{tin) is the initial value of the inflaton, and e is the average value of e in 
the slow roll regime. However, in this approximation, any deviation of e from a constant 
value is integrated over time, which can lead to signiflcant deviations in the evolution 
of (p towards the end of inflation. 

2.3. Linear perturbation equations 

The EM fleld does not contribute to the background expansion but comes into play at 
the perturbative level. We study the generation of perturbations in the FL background 
during inflation. We work in longitudinal gauge where the metric perturbation is 

dg^Jx^'dx" = a^{t){-2^dt'^ -2(!>6ijdx'dx^) (15) 

and the two scalar degrees of freedom, "$ and $, coincide with the gauge- invariant 
Bardeen potentials [13]. 

We work in Coulomb gauge throughout, i.e. A^ = (0, A^) with diA^ = 0. To lowest 
ordeil§|, the inhomogeneous Maxwell equation (JTj) with the axial current becomes 

i,-V% = -f{^)^ei,kdjAk. (16) 

where the Euclidean Laplacian is deflned as V^ = 6^^didj and the totally anti-symmetric 
symbol eijk satisfles 6123 = 1. Note that these equations are like in Minkowski space, 
there is no coupling to the scale factor. For a constant axial coupling, /'(</?) = 0, the 
sourced Maxwell equation reduces to the standard free wave equation and no fluctuations 
are amplifled during inflation. 

§ Since the EM energy density is quadratic in the fields, one considers the vector potential and the 
electric and magnetic fields to be at half order in linear perturbation theory. 



The evolution of perturbations in the scalar field is also altered by the axial coupling. 
At linear order, the scalar field equation (|6]) acquires a source term 

6(t) + 2H6(j) - V^6(j) + aV"((^)50 - v?(3<i> + ^) + 2aV'((^)^ 

= -a-^f'{ip)eijkAidjAk. (17) 

For any scenario where EM perturbations are generated, it is important to investigate 
the effect of this source term on the generation of scalar perturbations and through 
these on the primordial curvature perturbations. For instance, this has been studied in 
the case of natural infiation [16] . We leave a general discussion of such effects for future 
work [47] and concentrate here on the generation of magnetic fields. 



2.4- Physical properties of electromagnetism in the expanding universe 

The four-vector potential is generally covariant and its evolution is independent of 
the choice of coordinates. However, for an observer, the physical EM field manifests 
itself in terms of electric and magnetic fields which are intrinsically frame dependent 
quantities. Measured by an observer with four- velocity u^, with u'^Ua = —1, the electric 
and magnetic fields can be covariantly defined as 



E^ = F^^u'^ (18) 

B^ = ^v,.a^-yF'"'u'^ = F^aW". (19) 

These are both three- vector fields in the sense that they are orthogonal to the observer 
velocity, Eau"' = = Bau"'- In a perturbed FL metric an observer has the four- velocity 
u^ = a~^{l,0) + 0{1). As a consequence one finds (in Coulomb gauge) 

(E^) = (O, -^i.) , (B^) = (o, KjkdjAf}j (20) 

up to first order. With respect to an orthonormal basis comoving with the observer (or 
generally the Hubble fiow) , we define the Euclidean three- vector fields E and B through 

(E^) = afo, e), {B,) = afo, b) . (21) 



As discussed in Ref. [28], in a highly conducting plasma, magnetic fields should scale as 
a~^ with the expansion, while the electric field is damped away. The three-vector fields 
E and B defined above show exactly this property. We therefore rescale the fields by 
a factor of a^ such that the effect of the expansion is absorbed, i.e. B = a^B. In terms 
of the vector potential, the rescaled fields become 

Ei = —Ai , Bi = eijk djAk . (22) 

Using these expressions and the field equations for A^, one can derive Maxwell's 
equations for the rescaled fields E and B, which take the same form as in a Minkowski 
space-time. 

The physical properties of the EM fields can now be described in terms of the 
rescaled fields. Here, we are mainly interested in the energy density and the helicity 
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density of the fields generated during inflation. The energy density of the magnetic field 
is 

PB = \b^B^ = ^B B = ^a-^B-B = a-^pB (23) 

and analogously for the electric energy density. The helicity density is given as 

S) = AaB" = a~^5'^AiBj = a~^^. (24) 

(Deliberately, we are not denoting the vector-potential in bold face because Ai are the 
spatial components of the four-co- vector A^ in the covariant representation with respect 
to the coordinate basis {dx'^), while e.g. B are the components of the four-co- vector B 
with respect to the orthonormal basis {adx^) as discussed, for instance, in Ref. |28].) 

3. Electromagnetic quantum fluctuations 

To investigate the generation of EM fields during inflation, we consider the evolution of 
quantum fluctuations of the EM vector-potential. The coupling of the vector-potential to 
the background evolution of the inflaton via the inhomogeneous Maxwell equation (TT6|l 
can lead to the amplification of EM quantum fluctuations. The amplification depends on 
the axial coupling /(</?) and therefore on the evolution of ip. In Appendix A we review 
the quantisation of the vector-potential in an expanding background. The result is that 
Maxwell's equations lead to an evolution equation for the Fourier modes of the quantised 
field which reduces to the free wave equation in the absence of the axial coupling. We 
first discuss this evolution equation and then summarise the physical observables of the 
EM field expressed in terms of the solutions to the mode equations. 

3.1. Evolution of the Fourier modes of the vector potential 
We introduce the orthonormal spatial basis as 

(£'^,e^,k) with \e'[\^ = l,k = k/k, (25) 

and 

4 = i= {el ± tel) . (26) 

In radiation gauge, the vector potential then takes the form 

A = AiEi + ^2^2 = A+£+ + A-£- . (27) 

After quantisation of the vector potential, we can study the evolution of the Fourier 
modes, Ah(t, k), with respect to the hehcity basis for the polarisation states, h = ±. We 
find that the helicity modes satisfy the wave equation with a time dependent mass term 
corresponding to the modified Maxwell equation ( fT6|) for the classical vector-potential, 

A + [k^ + hkf{ip)ip] A = . (28) 

The fact that the sign of the respective helicity mode, /i = ±, appears explicitly in the 
evolution equation of the Fourier modes, leads to a different evolution of the two helicity 



9 

states and therefore, to the generation of magnetic hehcity. Also note that the scalar 
field couples to kAh and the coupling function itself, f'{ip)i^, only depends on time, as 
If is the background value of the inflaton. The solutions to this mode equation for a 
given coupling /(v?) fully determines the spectrum of the generated EM fields. 

Let us compare the mode equation ( l28l) to the non-helical case with the coupling 
f{ip)F'^. One obtains a similar evolution equation [28] 

A + 2^A + fc'A = 0. (29) 

Redefining the mode functions as Ah = fAh, one can rewrite this in the form 



Ah 



f 



Ah = 0. (30) 



We observe two significant differences to the helical case: firstly, the two helicity states 
couple with the same sign and therefore, no helicity is generated. Secondly, the scalar 
field couples to Ah (or Ah) as opposed to kAh in the helical case. As we shall see, this 
leads to significant differences in the spectra obtained for the two types of couplings. 
The reason is that the additional factor k in the helical case leads to a suppression of 
the coupling term at super-Hubble scales. 

To visualise this, let us compare the importance of the different terms in the helical 
mode equation fl28|) . We define Jn to be the logarithmic derivative of the coupling 
function, /, with respect to the scale factor 

f.^JL^iMp^ ±/'(^)^™p. (31) 

dm a H 

where A^ stands for the number of e-foldings and is defined as A^ = ln(a/ain). Note that 
/tv is the dimensionless part of the coupling term appearing in the mode equation ( 128|) . 
First consider super-Hubble scales, k <^'H, where we can approximate |c?t| ~ "H so that 
f l28|l reduces to 



1 + :^ + hfNj^ 



Ah ^ 0. (32) 



The second term is small by definition and the coupling term can be important on 
super-Hubble scales only if / is a rapidly varying function of time, i.e. /n = f /'H ^ 1 
and then only for as long as Jn > 'H/k. On sub-Hubble scales, for /c ^ "H, we may 
approximate Ah ~ k'^Ah- As long as /c :^ In'H the first and second terms in Eq. (128!) 
dominate. Hence the coupling term is typically relevant only at Hubble crossing k ^ T-L 
during a few Hubble times at best. In Fig. [H we show the evolution of the two helicity 
modes for a given wavenumber in case of a power law coupling, / oc ip^. Clearly, the 
modes only feel the axial coupling around horizon crossing. This fact turns out to be 
relevant for the resulting spectrum. We will discuss this issue in more detail in Sec. 14.21 



10 








1 1 1 1 1 


1 1 1 1 1 1 1 




4.6 






/-^^ 


- 


4.5 


- 






- 


2 4.4 


y^ 






- 


1 


\ 






- 


4.3 


N 


\ 
\ 




- 






\ 




- 






\ 
\ 




- 






\ 




- 


4.2 


- 


\ 




- 






N. 




- 








"-^^ 


" 






1 


1 1 1 1 


, ~i ~ r 



10 



12 

N 



14 



16 



18 



Figure 1. The evolution of the two hehcity modes (+ soHd red, — dashed blue) 
for wavenumber fc = 10/Mpc is shown as a function of the number of c-foldings, iV, 
during inflation. Here we consider an axial coupling function f oc tp^, as discussed in 
detail in Sec. 14.21 Both modes feel the axial coupling only around horizon crossing, at 
-Across — 10, while the evolution ceases and the modes saturate quickly after crossing. 
Here inflation ends at iV ~ 60. 



3.2. Vacuum solutions and initial conditions 

To study the generation of perturbations during inflation, we need to set tlie initial 
conditions when the mass term ± kTifN in Eq- ( l28l) can be neglected, i.e. well inside 
the Hubble horizon at early times. This condition is usually formulated in terms of 
the variable x = —kt which approaches infinity in this limit. During slow roll inflation 
H ~ — 1/t, so that X = —kt ~ k/Ti ^ 1 if the mode with wavenumber k is well inside 
the Hubble horizon. 

From the mode equation (l28l) . we see that if initially 

fwi'^m) = f'{^m)^m/'Hin < kfHin (33) 

the axial coupling term can be neglected with respect to the k"^ term and the mode 
equation becomes a free wave equation. Its solutions are plane waves. We match to the 



incoming vacuum solution described in Appendix A 

Ah{t,k) = A^''%t,k) = {2k)-^'^e-'^' 



for 



'kt > 1 



(34) 



This is used in the following as initial condition for the solutions of the full mode 
equation. Notice that the free plane wave solution only yields a valid initial condition 
if Eq. ( l33l) is satisfied. 
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3. 3. Power spectra and physical quantities 

The statistical distribution of the EM fields as seen by an observer can now be quantified 
in terms of a given solution for the helicity modes of the EM vector potential. We define 
the magnetic power spectrum and relate it to the magnetic energy and helicity density. 
If the magnetic field generated by some process is statistically homogeneous and 
isotropic, its spectrum is determined by two scalar functions Ps{k) and PA{k). Since 
the magnetic field is a divergence-free vector field the two-point function of the Fourier 
components of the magnetic field can be written as 

mt,k)B*{t,q)) = ^^6ik-q){iS,,-kk)Ps{t,k) (35) 

where Ps and Pa are the symmetric and anti-symmetric parts of the power spectrum, 
respectively. The symmetric part of the spectrum determines the energy density while 
the anti-symmetric part corresponds to the helicity density: 

{B,{t,k)B:{t,q)) = {27if6{k-q)Ps{t,k) (36) 

{Mt,k)B*{t,q)) = zk-^{kABUt,k)B*{t,q)) 

= k~\27if6{k-q)PA{t,k). (37) 

With respect to the helicity basis, see [Appendix A[ the spectra can directly be written 



as 

Ps/Ait,k) = kmA4t,k)\'±\A4t,k)\') , (38) 

where the upper sign corresponds to Ps and the lower sign to Pa- Here we use the non- 
trivial result widely applied in infiationary cosmology that at late times, the vacuum 
expectation values of the fields generated during infiation can be interpreted as stochastic 
power spectra. 

We define the magnetic energy density per logarithmic wave number via 

so that 

i|L(,,) . ^,Psit,t). (40) 

Similarly, we define the magnetic helicity per logarithmic wave number as i^ = A^B", 

it^k) =—PA{t,k). (41) 



d\nk^' ' 27r2' 

Finally, the electric field is given by the time derivative of the vector-potential and 
thus its contribution to the energy density, 'pe = EaE°'/2, is computed to be 

^(«.*)=(|ji(l^(*.*)P+l^-(a-)l^). (42) 

Notice that any electric fields produced during infiation will be damped very rapidly 
after infiation due to the huge conductivity of the primordial plasma. 
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4. Analytic solutions during slow roll inflation 

In this section, we first derive a condition on tlie coupling function f{ip) such that the 
theory can be treated perturbatively Even though we cannot prove that the theory 
does not make sense otherwise, our treatment is perturbative and so we can really 
trust it only if interactions are small. In our subsequent analysis in Section [6] we shall, 
however, analyse the results also if this condition is not satisfied. We then investigate 
two different forms of the coupling function and obtain upper bounds on its parameters 
without specifying a model of inflation. Finally, we solve the mode equation in slow roll 
for a constant Jn analytically, derive the resulting magnetic field power spectrum at the 
end of inflation and discuss the issue of backreaction on the background evolution of 
the infiaton. 

4.I. Condition on the coupling function 

In order to derive a condition on the coupling function such that the theory can be 
treated perturbatively, the interaction between the scalar field and the EM field should 
be small at all times. For this we demand that the ratio of the actions for the interaction 
term to the free EM term should be less than unity, i.e. 

Si[(p,A^] 



•Jcmi^/. 



< 1. (43) 

In [Appendix B| we show that in a FL background the free EM action can be written as 
Sem[A,] = V [dt{Mt)-Mt))v (44) 



where {■ ■ ■)v denotes the average over volume V in coordinate space and is assumed to 
be equivalent to the expectation value defined in Eq. fl39l) . Similarly, the interaction part 
of the action can be expressed in terms of the magnetic helicity as (see Appendix B[ ) 

1. 
2 



Si[(P,A,] = -V I dt f {S){t))v . (45) 



To compare the two parts of the action, one has to evaluate the time integrals for a 
given scenario. For the case of a maximally helical magnetic field, e.g. |^+| = |^| with 
|^_| = and assuming a power law form for the spectrum given by 

Ps = Pa = A; VI' = B^,k\ (46) 

we find that 

where k^a^{t) is the cut-off scale which is the smallest scale crossing the Hubble scale 
at time t when the volume average is computed, hence fcmax — 'H(t). To estimate the 
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energy density for the electric field, we can approximate |^| ~ "Hj^l at large scales 
which then leads to 

After evaluating the volume averages of the two energy densities and the helicity, we can 
now compare the two actions. Ignoring numerical factors, the comparison condition (l43ll 
can be written as 



Si[4>,A,] 



•Jcmi^^J 



Jdtn 



n+3 



< 1. (50) 



Now, for the above inequality to be satisfied, it is sufficient to require |/| < Ti or 
equivalently |/Ar| < 1. We believe that this condition is quite generic and does not 
depend on the functional form of the coupling function. The above analysis indicates 
the fact that the condition for the perturbative treatment of the theory does not depend 
on the coupling function, but on its derivative. This is not surprising as a constant / 
only yields a surface term which does not affect the dynamics. 

In fiat spacetime the above condition has to be replaced by |/| < /cmax, where /Cmax 
denotes the UV cutoff of the modes under consideration. 

4-2. Power law and exponential coupling 

We first study the special case of a power law coupling, 

fiv) = /o (^Y . (51) 

\mpj 

Here /o and p are constants. The sign of /o determines which helicity is amplified by the 
coupling and, thus, we can take /o to be positive without loss of generality. Let us first 
investigate bounds on the parameter space of /o and p under the condition \f^\ < 1. 

During slow roll inflation, we use Eq. (fT3l) to express the coupling term in the mode 
equation (l28l) . /at = f'{ip)ip/'H, in terms of the slow roll parameter e 



h ^ ±fopV2e{^\ (52) 

\mpj 

where, in principle, e is a slowly varying function of time. If we want to satisfy the 
condition of perturbativity, |/Ar| < 1, the above equation leads to an upper bound on 
the value of /o as 

/„</.- = A.»i4^(^y"|. (53) 



For a given inflationary scenario, one can invert the definition of e to express the above 
bound only as a function of (p. Though it turns out to be easy for large field inflationary 
models, it is not so straightforward in the case of typical models of small field inflation. 
Once the upper bound on /o is known, one can calculate the maximal possible value of 
the coupling term i.e. /^^^(p) = /Ar(/™'^''(p),p)- 
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As a second case, we consider the example of an exponential cliiral coupling of the 
following form 

/(v.) = /oe'^^/'"- (54) 

where /o > and a are constants. Analogously to the power law case above, we calculate 
the coupling term /tv in terms of e which is given by 

JN ^ ±/o«y2ie-^/'"^ (55) 

We again apply the constraint |/Ar| < 1 to find the bound on /q: 

/o < fr"" = ^min(^e-"^/-4. (56) 

As before, once Z™'*^ is known, the maximum value of the coupling term f™'"^^ can be 
calculated. 

In both models, the power law and the exponential chiral coupling, the constraint 
I/tvI < 1 leads to an upper bound on the overall amplitude of the coupling. The time 
dependence of /at, on the other hand, depends on the details of the inflation model. 
Though, generally we remark that for a given mode k the coupling is only active for 
a few e-foldings during which the mode crosses the Hubble scale, and moreover, if f^ 
had a significant time variation its overall amplitude would be suppressed very strongly 
by the bound on /o < f^^^. Therefore, we conclude that /at can safely be considered 
roughly constant (at least during horizon crossing). 

4-3. Analytic solution for a constant coupling term 

We now solve the mode equation fl28l) for the case where f^ is approximately constant. 
The solution presented here will be valid for any constant value of Jn- In the slow roll 
regime, the full coupling term is, f'{(p)'f = /n'H ~ —fN/t, i.e. inversely proportional to 
conformal time and, consequently, the mode equation can be solved analytically. It is 
convenient to use x = —kt as the time evolution variable. For each scale, k, the initial 
condition in the asymptotic past is set well inside the horizon, i.e. for x — ?■ oo, while 
inflation is considered to end when x — )■ 0. With a prime denoting the derivative with 
respect to x, the mode equation reads 

A'lix) + [1 + hf^x-') Ahix) = 0. (57) 

The free solution is a good approximation at early times, x ^ IInI- The general solution 
to the mode equation (1571) is [19] 

A(x) = C^Goiy,x) + C2Fo{y,x) (58) 

where y = —hfN/2 and Go and Fq are the irregular and regular Coulomb wave functions 
of order zero, respectively. As initial condition we require the solution to approach the 
free solution, Eq. (!34l) . for x — ;■ oo. It turns out that the combination Gq ± iFo has the 
desired limit, as given in Ref. [49] : 

Go{y,x)±iFo{y,x) — > exp±i[x - y\n{2x) + ao{y)] (59) 
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with cro(y) = argr(l + iy) = jeV + 0{y^) and 'Je — 0.58 being the Euler's constant. 
For large x, the second term in the exponential, y\n{2x), can be neglected with respect 
to X. Comparison with the free solution shows that the plus sign corresponds to the 
incoming vacuum solution and we have to choose the initial amplitude 

^(Ci-^Ca) = (2fc)-^/2 ^^^ ^(^Ci + tC2) = 0. (60) 

The factor with ctq only acts as an overall phase and has no physical significance. The 
normalised full solutions are 

Ah{x) = {2k)-'/^[Go{-hh/2,x)+tFo{-hfN/2,x)] (61) 



for h = ±1. This solution of the mode equation has also been found in Ref. [32] for 
A^-flation. Below we shall argue that it is very general. 

To understand the effect of the axial coupling on the growth of EM quantum 
fluctuations and to compute the magnetic power spectra at the end of inflation, we 
analyse the late time limit of this solution, i.e. when x — )■ 0. Using the approximate 



expressions derived in Appendix C , we find the asymptotic limit at late times as 



A(x) -^ i2k) 



-1/2 



exp(-/i7r/Ar/2) sinh(7r/7v/2) 1 ^^^ 

^n^ • ^^^' 

t^Jn/2 

The late time behaviour of both helicity modes is independent of x, and therefore of t, 
and the scale-dependence is not changed with respect to the free solution. The modes 
are coherently amplified while crossing the horizon, before they saturate outside the 
horizon. Notice also that we did not use the asymptotic limit given in 14.6.9 — 10 of 
Ref. |l9] because it is not correct. (In this context, see discussion in Appendix C[ ) 



Given the late time behaviour of the helicity mode functions, it is easy to compute 
the magnetic field power spectra produced at the end of inflation. The symmetric power 
spectrum is 

while the antisymmetric one is 

PA{k) = fcS^!M![M^ (64) 

Notice that both spectra are proportional to k and only their amplitude changes with 
the coupling strength /^y. In Fig. [21 we illustrate the /c-independent amplification factors 
of Ps and Pa as a function of /at. The larger |/Ar| the smaller the difference between Ps 
and Pai i.e. the more helical the magnetic fields become. As one infers from Eqs. (!63|l 
and ( 16^ . both amplification factors tend to {2nfN)~^ exp(7r/Ar) for large values of /at. 
The magnetic energy density per logarithmic wave number at the end of inflation 
can directly be computed to be 

dpB ,. ,^ k^ sinh(7r/jv) _ k^ ^ 

-{tcnd,k} = — — — = -^—S [Jn). (65) 



dink ' a^j^d An^f, 



N "g^j 
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Figure 2. The symmetric (solid, red) and antisymmetric (dashed, blue) power spectra 
of the magnetic fields, Ps^A/k, are plotted as a function of the effective coupling 
constant |/Ar|. For vanishing coupling, the vacuum solution, Ps/k = 1 and PA/k = 
is recovered, while the larger the coupling, the smaller the difference between Ps and 
Pa- As discussed in the main text, for the theory to be perturbative we must require 
I /at I < 1 and therefore, the amplification of the magnetic fields is small. 



Here, we define 5^ (/at) to be the amplitude of magnetic energy density spectrum. 
Since its spectrum is blue, the magnetic energy density is dominated by the cut-off 
scale, which is set by the last scale that exits the horizon before the end of inflation, 
kc = "Hend = {ciH)end- With this we obtain 

PB(te„d) ^ \h^^^S\M. (66) 

By means of the Friedmann equation, we have 



i^B tend = ^ 10 • ^" 

Ptot 12 \mp J 

The generic condition for backreaction to be negligible then is 

S{fN) < ^. (68) 

-"end 

If we assume that the reheating phase is short so that SrupH^^^ = ptot — T^, where T* 
is the reheating temperature, we can write the above bound from backreaction as 

2 

SUn) < ^. (69) 

For couplings that are perturbative in the naive sense discussed above, |/Ar| < 1, S^Jn) 
is of order unity and, therefore, no backreaction on the background evolution is expected 
if inflation ends well below the Planck scale. However, the results above hold for any 
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constant value of /at, even if it was larger than unity. Since we cannot prove that 
perturbativity is strictly required, we shall therefore not restrict ourselves to this case 
in the following. For instance, if the EM field was coupled to a large number, A^, of 
pseudo-scalar fields then /tv — ;■ \/JJ Jn, a situation where the effective coupling term 
can be large without spoiling the perturbativity of the interactions with the individual 
scalar fields, see Ref. |42] . 

In the more general case where /tv is time dependent, in principle we have to 
evaluate f^ at horizon crossing. This would lead to a slight modification of the spectrum 
but would not spoil the discussion of backreaction above. In the following we neglect 
this effect, as it is quite irrelevant for the few orders of magnitude in k which we are 
interested in, see Fig. |H 

At the end of inflation and after reheating, we expect the universe to be filled with 
relativistic standard model particles, a relativistic highly conducting plasma. In this 
medium, the MHD approximation is valid and electric fields are rapidly damped away. 
We therefore do not discuss the electric field spectrum which will not survive reheating. 

The helicity density per logarithmic wave number at the end of inflation is simply 
given by 

dS) k^ cosh(7r/^) - 1 

dink a^^d 27rV7v 

Note that S^ oc ps/kc, which therefore has to be constant when helicity is conserved. 

We believe that this result is more general than the particular cases studied here: 
whenever the function /^r is rapidly varying so that /Ar(t) = constant during slow roll 
is no longer a good approximation, the fact that we require |/Ar| < 1 during the entire 
period of inflation means that f^ must be oscillating. Though, because the coupling is 
active only for the small number of e-foldings during which a mode crosses the Hubble 
scale, resonant amplification of modes at horizon-crossing is not likely to be efficient. 
We have checked this statement numerically using different forms of oscillating coupling 
terms. We, therefore, conclude that as long as the slow roll approximation is valid and 
the theory can be treated perturbatively, the amplification of helical magnetic fields is 
always mode independent and consequently leads to a n = 1 spectrum for both the 
magnetic field and the helicity. 

Note that this result differs significantly from the non-helical case. There, the source 
term in the mode equation is of the form /// which is typically oc 1/t^. The solutions 
are then Bessel functions and the Bessel function index, which determines the spectral 
index at late times, is related to the (nearly arbitrary) pre-factor. The difference comes 
from the fact that a term of order 1/t^ is relevant during all the time when the mode is 
super-Hubble, —kt < 1, while a term of order {k/t)Ah is relevant only around horizon 
crossing. On super-Hubble scales, it is dominated by the Ah ~ Ah/t"^ term while on 
sub- Hubble scales, the k'^Ah term becomes dominant. 
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5. Deviations from slow roll 

In the previous section, we discuss two different functional forms of the axial coupling of 
the infiaton to the EM field, namely a power law and an exponential. Within the slow 
roll approximation we find that Jn can safely be considered constant. This always leads 
to a magnetic field power spectrum proportional to k. In this section, we first confirm 
our analytical findings numerically. Second, we explore the possibility of obtaining 
a different magnetic field spectrum by introducing a short deviation from slow roll, 
motivated by the fact that such deviations can provide a considerably better fit to the 
angular power spectrum of the CMB anisotropies than the predictions from typical 
single field inflation models, see for instance Refs. [501 EB Ell [531 EH 155] . 

To compare the analytical result, Eq. ( 163|) . to a full numerical solution, we solve the 
background evolution of ip with a quadratic potential, Vo{(p) = \rn?ip'^, and integrate 
the evolution of the modes Ah{t, k) to compute the magnetic power spectrum Ps{k) at 
the end of inflation. 

A short deviation from slow roll can, for example, be achieved by introducing a 
step in the quadratic infiaton potential as follows [50l EH E2l l56] 

'ip-ifQ^ 



y^i^) = 2"^v^ 



/3tanh 



(71) 



Av9 

Here /3, ipQ and A^p characterise the height, the location and the width of the step, 
respectively. Such a deviation from slow roll, in general, leads to a burst of oscillations 
in the primordial power spectrum of curvature perturbations. In what follows, we turn 
our attention to the possible effects of such a deviation from slow roll on the magnetic 
field power spectrum. 

We perform the comparison with a power law coupling function, /(v?) = fo{ip/mpy. 
In this case we find an exact expression for the slow roll parameter: y/2e = 2'mp/ip. 
Using this in Eq. ( 153|) the upper bound on the parameter /q becomes 

^ J_( {Vrajrupf-^ for p > 2 

° 2b| \ {^^Jmpf-^ for p < 2 

and v^min <'~p< V^max- With /o fixed to be f^''''{p) the coupling term /^^^ = fNif^^'^'^p) 
reads 

f (<^/<^max)^"^ for |) > 2 
/r" = < (73) 

[ (<^/<^min)^"^ f0Tp<2 

Thus, the most simple choice of parameters is: p = 2 such that f^^^ = 1/4 and f^^^ 
is exactly constant in slow roll. Furthermore, we choose /o = (4/5) Z™*^^ = 1/5 to make 
sure that also in the case with the deviation from slow roll the condition |/jv| < 1 is 
respected. 

As we discuss earlier, the coupling term f^ is typically relevant only around Hubble 
crossing for a few Hubble times. The coupling term is, in general, a function of time but 
this time dependence can be translated into a scale dependence by identifying a time 



^ 
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Figure 3. The coupling term I/atJ at Hubble crossing of the wave number k is plotted 
for the power law coupling with p = 2, fo = 1/5. The horizontal dotted blue line 
represents the constant value of /jv = 4/5 in the slow roll approximation. The dashed 
black line indicates the behaviour of the coupling term for the potential in Eq. ([7T|) 
with /3 = 0, while the solid red line is the same for (3^0. For the later case, we have 
used the best fit values of the parameters of the potential as in Ref. [55] . The bump in 
the coupling term for /3 7^ arises due to the short period of deviation from slow roll. 



with the corresponding Hubble crossing scale, k = 'H{t). In Fig. |3]we plot the coupling 
term /^r at Hubble crossing of the mode k for the power law coupling with p = 2 and 
/o = 1/5. Comparing the slow roll value, Jn = 4/5, with the numerical results for the 
slowly rolling model and the case departing from slow roll, it is evident from the figure 
that a deviation from slow roll leads to a bump in the coupling term as compared to 
the slow roll case and therefore, one can expect an effect in the magnetic field power 
spectrum on the scales which exit the Hubble radius around the time when the bump 
in the coupling term occurs. 

In Fig. m we plot the relative deviation of the magnetic field power spectrum over 
an exact /c-spectrum as a function of k for the power law coupling. The exact numerical 
solution for slow roll deviates slightly from the Ps oc k spectrum due to the slight scale 
dependence of iS^(/Ar). Modulations in the spectrum arise as a result of a deviation from 
slow roll. We find that for the best fit values of the parameters of the potential (17T]1 . 
the spectrum of the magnetic field is not strongly modified. Indeed, we conclude that 
even the deviation from slow roll, within the limits required by CMB data, does not 
significantly modify the magnetic field spectrum. We find that the exponential coupling 
leads to a similar behaviour for the coupling term and the magnetic field spectrum. 
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Figure 4. The relative deviation of the magnetic field power spectrum from an exact 
fc-spectrum is plotted as a function of the wave number k for the power law coupling 
with p = 2, /o = l/5. The dashed black line indicates the numerical solution for slow 
roll inflation (/3 = 0) while the solid red line is the spectrum resulting from a deviation 
from slow roll. The horizontal dotted blue line is the magnetic field spectrum from the 
slow roll approximation, /jv = 4/5. The deviation from this approximated spectrum 
is always small. Only the scales which exit the Hubble radius around the time when 
a deviation from slow roll or equivalently a bump in the coupling term occurs are 
affected more strongly. 



6. The magnetic field at the end of inflation and its further evolution 

Our main results are Eqs. (pSj) and f l70p which determine the magnetic energy density 
and helicity density at the end of inflation. Note, however, that we do not renormalise 
the energy density. Hence even for /at = 0, we obtain the non- vanishing result 



dpB 



(^end; "^j 



k' 



which comes purely from (not amplifled) vacuum fluctuations and may be considered 
unphysical. However, from Eq. f p8|) it is clear that Ps{k) > |Pa(^)| by deflnition. 
Hence the physical result cannot be obtained by a simple subtraction of the vacuum 
contribution as then Ps oc /at^ would become smaller than P4 oc f^ for small values 
of /tv, see Eqs. ( 163|) and ( 164|) . On the other hand, for /jv > 1, the vacuum contribution 
becomes subdominant and it is no longer important to subtract it. We shall therefore 
not perform any renormalisation of the magnetic energy density but just keep in mind 
that our result becomes dominated by vacuum fluctuations in the limit /at — )> 0. 
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After inflation, the thermal cosmic plasma contains many relativistic charged 
particles and can be treated as an MHD plasma. During the process of reheating, 
the Reynolds number becomes very high and MHD turbulence develops. In the MHD 
limit the electric fleld is damped away and the magnetic fleld evolves by two different 
processes: it is damped on small scales and it undergoes an inverse cascade due to 
helicity conservation [31]. Numerical studies have shown that very soon damping on 
small scales leads to a maximally helical field (for which either Aj^. or ^„ vanishes) which 
then continues to evolve via an inverse cascade, see [32l|33]. The inverse cascade is active 
as long as the Reynolds number of the cosmic fluid at the scale under consideration is 
larger than one and the fluid is therefore turbulent [57]. The damping scale fcdiss(^) is 
the scale at which the Reynolds number becomes of order unity. On scales smaller than 
^diss(^) the magnetic fleld and the turbulent motion of the fluid are damped exponentially 
by viscosity. 

In the following we investigate how the spectrum of helical magnetic flelds evolves 
during the turbulent epoch. We flrst discuss the evolution of the correlation scale of the 
magnetic fleld and the duration of the turbulent phase, before computing the magnetic 
energy spectrum at the end of the inverse cascade. Here we assume that the reheating 
epoch is relatively short and ends at t*. This corresponds to the reheating temperature 
T*, and because of radiation domination, we can approximately use t/t^, ^ T^^/T during 
the turbulent phase. 

The helical magnetic fleld from inflation always has a blue spectrum and is therefore 
dominated by the largest wavenumber crossing the Hubble scale at the end of inflation 
or, for simplicity, reheating 

psiQ ^ k^S'^ifN), with k^^n^. 

Accordingly, the correlation scale (which is roughly given by the scale at which the power 
spectrum peaks) is initially kc{t^) = k^. As an example we consider T^ = lO^^GeV and 
flnd 



f^cyt*) — f^* — rL* — ii^a,^ — \ n, 9o 9* -'* 

V 6 nip 



10-^^GeV( 



g* \ 1/6 f Z 



200/ VlOi^GeV 



.2xlO-Mpc-(|^)'"(^^). (74, 

Here asB is the Stefan-Boltzmann constant, asB = vr^/15 in our units, g^, and Qq = 2 
denote the number of relativistic degrees of freedom at t* and today, respectively, and 
To = 2.73 K is the present CMB temperature. 

Let us assume that the inverse cascade starts at t*. In Ref. [33] it was found that 
during the inverse cascade of a maximally helical magnetic fleld the total rescaled energy 
density scales like 

PB{t) oc (t/t*)-'/^ 
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Figure 5. The Reynolds number Re(r) (black solid) and the comoving correlation 
scale fcc(T)/fcfin (red dashed) are shown in log scale as a function of log(T*/r) with 
T* = 10^"' GeV. Turbulence and with it the inverse cascade terminate at Tfin — 1 GeV 
and fcfin = fcc(ifin) =i (10-12 Mpc)-i. 



and the comoving correlation scale evolves in the same way 

such that the ratio ps/kc which is proportional to the rescaled helicity density remains 
constant. This continues until tfin, the time when the damping scale has grown up to 
the correlation scale, kdissit&n) = ^c(^fin)- After tgn the inverse cascade and turbulence 
cease and the magnetic field evolves solely by fiux conservation on large scales 



dp 



B 



ik,t) 



dp 



B 



{k,t 



fin J 



for t > tfin , k< k^i^^it) 



dliik dliik 

and viscosity damping on small scales, k > kdiss{t). At the end of the inverse cascade, 
the correlation scale of magnetic field spectrum has moved to 

^fin = ^c(^fin) = ^* (^fin/^*) 

and the total energy density is reduced by the same factor. 

To compute the ratio Tfin/T* we need to determine the temperature at which the 
Reynolds number becomes unity. In Appendix A of Ref. [M] the Reynolds number at 
very high temperatures is estimated to be 



Re{k,T) oc 



aT lPB{k) 



k 



Pf 



for a given scale k. Here pf is the energy density of the fiuid which contributes to the 
turbulent motion. In perfect thermal equilibrium pf = p. More precisely, denoting the 
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Reynolds number at the beginning of the inverse cascade by i?* = Re( 


%{Q,T,), it is 


found 




r R4^Y^' for T>%^- 


= 100 GeV 


Re{K{T),T) = \ . xi/3. xii/3 

RJ^) (t-) for T < Te„ . 


(75) 



Setting g, ~ 200 and T, = 10^^ GeV yields R, - a,Tjk,{Q - ^^/{asBQ*) mp/T, 
O(IO^). (For simplicity we have set Pb/ Pf ~ 1 foi' this value.) Note that after inflation, 
until the electroweak transition at Tc„, the Reynolds number of the fluid is actually 
increasing. This comes from the fact that it is inversely proportional to the comoving 
mean free path which is constant at early times. After the electroweak phase transition, 
collisions are much more strongly suppressed and the comoving mean free path grows 
like a^, hence the Reynolds number decreases rapidly. For more details see Ref. [M]. In 
Fig. [5] we show the evolution of both, the Reynolds number and the correlation scale of 
the helical magnetic field through the inverse cascade. 

Finally, we can now derive the generic scaling of Tgn with the initial temperature 
T^. With the scaling /cgn = A;*(Tfin/T^,)^/^ and the help of Eq. fl75|) for the evolution of 
the Reynolds number, we find that the Reynolds number becomes unity and the inverse 
cascade stops at Tgn given by 

— ~ 10 ( 1 , for T, > Te^ (76) 



so that 



5/11 



■fi,, ^ 10^'Mpc-M— f^) . (77) 



1014 GeV 

For T* = 10^4 GeV we obtain Tgn — 1 GeV and the correlation scale moves by about 9 
orders of magnitude from k^, to k^^, see also Fig. |5l 

Now we can trace the magnetic field spectrum through the inverse cascade. The 
spectral shape on large scales remains unchanged [33] • At late time, t > tgn we therefore 
obtain 

^f^u M ^ / ^(^-M ff)'te)' for k<k,^ 



dink 




else. 



With T*/Tfin ~ 10^4 this yields an amplification of the initial amplitude by 28 orders of 
magnitude on scales larger than l/Zcgn; see the sketch in Fig. O 

For a generic spectral index n, Eq. ( 178|) is replaced by, see [31] , 

{/ \ n+3 / \ 2(n+2)/3 
^(t*, k,) (^A j (^^J for k<k^^,t> tfi, ^^g^ 

else. 

Clearly, the smaller n the less significant is the amplification by the inverse cascade and 
for n = —2 there is no amplification at all. For n < —2 the above result does not apply, 
see [31. 
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In our case, where n = 1, Eq. (165 p can be written as 



-iU^K) ^ 4 X 10~''GeY'S\M [^f j^- 



rfln/t' '^' *' '■"'' V200/ VlO^^GeV 

with which we arrive at 

where we assumed g* ~ 200. This determines the final strength of the magnetic field on 
large scaleqJlJ, 

BW.3xlO-Gauss5(/.)(-) (^^.j^) (80) 

for k < fcfin. With the help of Eq. ( 1771) this can also be written as 



B{k) ~ 3x 10"^^ Gauss 5(/, 



N) 



k W T. ^'/'' 



lOiVMpcy VlO^^GeV 

for k < kfin- 

After the end of the turbulent phase, magnetic fields are damped on small scales 
by viscosity and evolve by fiux conservation, so that B = const, on large scales. For 
our typical value of T^, ~ lO^^GeV hence k^^ ~ 10^^/Mpc, for cosmologically interesting 
scales of A; ~ 10/Mpc the magnetic field is of the order of B{k = 10/Mpc) ~ 10"^" Gauss. 
This is much too small for dynamo amplification. For smaller reheating temperatures, 
T*, the Reynolds number grows less strongly and turbulence and the associated inverse 
cascade are of shorter duration. Therefore the value of B{k) at fixed k < /cgn is 
actually smaller for T* < 10^"^ GeV even though A;* is larger for a smaller reheating 
temperature, see Eq. (IHT]) . Considering the lowest value for which our treatment is 
valid, T, ~ Te„ ^ 100 GeV we arrive at Tgn ~ 8MeV and A:fin ^ 4 x lOVMpc but the 
magnetic field is only 

~ f k\^ 

B{k) ~ 6 X 10"^° Gauss ^(/tv) (-;— ) , for A; < kf^^ ~ lOVMpc . (82) 

At scales of 0.1 Mpc, this field is by far insufficient for subsequent dynamo 
amplification which requires seed fields of the order of at least 10~^° Gauss [31]. For 
an arbitrary reheating temperature T^ we obtain from (ISTl) 

B{k = 10/Mpc) ^ 3 X 10-^1 Gauss S{M f j^n^ 1 (83) 

Let us investigate how large this field can become in the best case in which our treatment 
may apply. Certainly we want to require backreaction to be unimportant, but we do 
not insist that f^ be small. Using (169|1 this requires 5(/Ar) < {mp/T^Y and hence 

S(A; = 10/Mpc) < 10"^' Gauss ( J (84) 

II The magnetic field strength is B = y^SirpB in Gauss units, where GeV ~ 1.4 x 10"'^^ Gauss, while in 
Heaviside-Lorentz units we have B = \J2pB and GbV^/a/Itt ~ 1.4 x 10^^ Gauss. 



25 



•30 - 



-40 



Of) 

o 



^ -50 - 



13 



-60 



-70 - 




Figure 6. The process of inverse cascade is depicted. We show ■j\^{kcit),t) as a 
function of kc{t) (dashed hne) as well as ^^(fc,^*), ^^ik,ti) and ^^:(fc,ifin) as 
functions of k (thick solid lines). We choose ti/t* = 10^. The parts of the curves to the 
right of the maximum are not reliable. The vertical line indicates the total amplification 
factor which is constant for fc < feg = /c(ifin)- The horizontal line indicates the amount 
by which the correlation scale increases during the inverse cascade. The energy density 
is in units of Gauss for S{fN) = 1 and k is in units of Mpc^ . 



To achieve a minimal necessary field for dynamo amplification of i? ~ 10"^*^ Gauss we 
would need T^ < 10^ GeV, a somewhat low inflation scale, but not excluded. 

Until recombination, fields on scales smaller than about 0.1(i?/10~^Gauss) Mpc 
are damped by viscosity [58]. The scales A; > 10(10~^Gauss/i?)/Mpc therefore do not 
survive the linear regime and will not be amplified before being damped. But even on 
these smallest "surviving scales" the magnetic field generated is in most cases too weak 
for dynamo amplification. 



7. Conclusions 

In this work, we have studied the generation of helical magnetic fields during inflation 
by adding a parity- violating term of the form f{(f))FF to the action. For two specific 
choices of the coupling function f{(f)), namely a power law and an exponential form, we 
found that, during slow roll inflation, the power spectrum of the generated magnetic 
fields is always blue with spectral index n = 1. The aspect responsible for this result 
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is that the hehcal couphng term in the evolution equation for the Fourier modes of the 
vector potential is active only around horizon crossing of a given mode and it is not 
varying rapidly during slow roll inflation. As a consequence, in slow roll inflation all 
Fourier modes are amplified in the same way and the final shape of the power spectrum 
is unchanged with respect to the initial one. 

We have also studied the effects of a short deviation from slow roll. If kept within 
the bounds allowed by the CMB data, such a departure from slow roll does not strongly 
modify the overall shape of the magnetic field power spectrum. In principle, a very steep 
coupling function, e.g. a double-exponential /(0) = /oexp[exp(a0/mp)], can lead to a 
different spectrum, n ^ \. However, if we want to satisfy perturbativity, |/Ar| < 1, for 
such a behaviour to be valid throughout inflation, a ridiculously small pre-factor /o is 
required. In the absence of a convincing physical motivation, we have not investigated 
such an extreme case. 

After inflation and reheating, an inverse cascade sets in and moves power from 
the correlation scale fc* ~ H* to larger scales, which is schematically represented in 
Fig. ini We find that typical values for the reheating temperature, e.g. T* ~ 10^^ GeV, 
can only lead to magnetic fields of the order of 10"^*^ Gauss on scales of 0.1 Mpc and 
smaller values are obtained for smaller reheating temperatures. Unless S{fi\i) ^ 1 and 
the reheating temperature is sufficiently low, T^, < lO^GeV, these field amplitudes are 
largely insufficient for dynamo amplification |31j . 

This result is quite generic. A power law or an exponential coupling to the FF 
term during slow roll infiation will generate sufficiently strong seeds for the large scale 
magnetic fields in galaxies and clusters without also leading to backreaction only if the 
infiation scale is sufficiently low, of the order of the electroweak scale. Even though for 
helical fields an inverse cascade does move power to larger scales, this is still insufficient 
in most cases. For typical infiation scales, to obtain sufficient fields on large scales after 
the inverse cascade, they must have been so large on small scales after infiation that 
backreaction cannot be neglected. Another possibility might be a very steep coupling 
which can lead to a different, n ^ 1 spectrum. However, only if the spectrum is close 
to scale invariant, it can have significant amplitudes on large scales without significant 
backreaction from small scales. 

Scenarios of infiationary magnetogenesis are often constrained by requiring that 
the backreaction of the generated magnetic field on the background evolution is small. 
Since the perturbations in the scalar field are affected by the presence of a non-minimal 
coupling as indicated in Eq. f lT7|) . it will be interesting to study the backreaction effects 
of the EM perturbations on the evolution of the primordial curvature perturbations. 
This may provide another tool to constrain inflationary scenarios of magnetogenesis. 
This issue will be addressed in a future project 
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Appendix A. Quantisation of the vector potential 

Always working in Coulomb gauge, we promote the vector potential Ai to a quantum 
mechanical operator, define conjugate momentum as W = 6S/6Ai, and impose the 
commutation relation 

[A{t,x), W{t,y)] = i5i^{x-y). (A.l) 

Here the transversal Dirac delta, 5;^^, is defined as 

with ki = ki/k and k = \k\. Let us compute the conjugate momentum in terms of 
the vector potential which we expand in terms of creation and annihilation operators 
in Fourier space. The canonical commutation relations of the creation and annihilation 
operators, together with the above commutation relation of Ai and W will lead to 
a Wronskian normalisation condition on the mode functions of the quantised vector 
potential. 

The conjugate momenta in the perturbed FL metric are 

m = a^g'^ (i, + f{ip)eJ^n^„^A^) . (A.3) 

The non-standard second term arises due to the axial coupling of the EM field to the 
infiaton. Since it involves the curl of the vector potential it is convenient to expand 
the operators in Fourier space with respect to a helicity basis. For each comoving wave 
vector. A;, we define the comoving right-handed orthonormal basis 



with 



and 



el =a(l,0), e'fik) ^ a(o, 6^2) , ej(fc) = a(o, fc) (A.4) 

A;.£^2 = = £^s^ e^Ae^ = k. lef |' = 1 (A.5) 



e^" = -et, €2-^ = 4- (A.6) 

The two transverse directions are combined to the helicity (or circular) directions 

4 = ^ (4 ± ^4) • (A.7) 
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(A.8) 

(A.9) 

(A.IO) 

(A.ll) 

_ . ,_ (A-12) 

h 

where h,h' G {+, — } and the star denotes complex conjugation. We now expand the 

vector potential in Fourier space in the helicity basis. 



These have the following 


useful properties 


k* k 






et = -e 


.k 

T 




ikAe'l = 


±4 




k _ k* _ 


Shh' 




\ ^ fc k 

2_^ ^h,i ^h,j 


= Sij 


rvirvj 



+ ef(fc)6t(fc)A*(t,A:)e-*-} (A.13) 

where the creation and annihilation operators, bl{k) and bh{k), satisfy the canonical 
commutation relations 

bh{k),bl{q)j = {27if5\k-q)5yM (A.14) 

Mk),h,.{q)] = = [biik),bl{q)] . (A.15) 

The vacuum, |0), is defined by bh{k)\0) =0. As a consequence of the commutation 
relation (lA.ip of Ai and its conjugate momentum, II-^, as given in (jA.Sp . the mode 
function must satisfy the Wronskian normalisation condition. In terms of the rescaled 
Fourier modes of the vector potential, Ah{t, k) = aAh{t, k), this reads 

Ahit,k)Alit,k)-AUt,k)Ah{t,k) = %. (A.16) 

Finally, the evolution equation for the Fourier modes can now be derived from the forced 
wave equation of the vector potential, Eq. (TT6l) . 

A(t, k) + {k^ + hkr{^)^\ A(t, k) = (A.17) 

where h = ± reflects the two helicity modes. 

Appendix B. Perturbative chiral interaction 

In order to derive a bound on the parameter space of the coupling function /(0), be it 
a power law, an exponential or any other functional form, we demand the interaction 
between the scalar field and electromagnetism to be small in order for perturbative 
quantum field theory to make sense. For example, we want the interaction picture to 
apply and the first terms in the Feynman-Dyson series for the scattering matrix to be 
small. Therefore we demand that |S'/[0, A^]| < 1 5*0111 [^/x]|- 

Let us express the free EM Lagrangian density in terms of the electric and magnetic 
energy densities (as defined in Sec. 12. 4p 

C.m = -\f' = 1{E'-B') = PE-PB (B.l) 
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and compute the free EM action in the FL background: 

5'em[-4^] = / d'^Xy/^Ccm = / d'^X O^ {pE - Pb) 

V / dt {Mt)-pB{i))v 



dt I d^x (pi 



PB) 



fB.2) 



where we rescaled the energy densities as ps = a^pB and defined the volume average 
(. . .)v in the last step. We assume that for quantum fluctuations the volume average 
equals the vacuum expectation value, (...), such that 



{Mt)h 



(Mt)) 



V 






~ k^dk 

^ k^dk 



' A+{t,k)^ + \A^{t,k)\' 


(B.3) 


'\A+{t,k)\'+ A.{t,k)^ 


, (B.4) 



'0 (27r)^ 
as we show in Sec. 13. 3[ 

To be able to compute the contribution of the interaction f{(f))FF to the action 
we pass to the language of differential forms by writing the electromagnetic 2-form 
F = ^F^^dx^'dx" = dA = di^A^dx^"), such that 

FFy/^d^x = 2FAF. (B.5) 

Then we use the facts that F = dA and dF = d{dA) = to integrate by parts 



SilAA 



5 ^-"MJ 



^ / d'x 



- I f{<P)dAAF 



gfi<t>)FF 

= ^j (d[fAAF]-dfAAAF^ = ^ j df A F A A . (B.6) 

In the last step we assume the boundary terms to vanish. Clearly, S'/[0, A^] vanishes for 
a constant /(0). Next, we express F in terms of E and B as measured by an observer 
with timelike four-velocity u, see Sec. 12. 4[ 

F = uAE-2(^^AB). (B.7) 

Using this we find 

dfAFAA = df AuAEAA + 2f^^u^''Bl^^Ap^^d^x. (B.8) 

In a FL background, the only valid observer is the Hubble fiow, u = —adt, and / is only 
a function of time, i.e. df = fdt, such that df Au = 0, f^au" = a~^f and f^aB^ = 0. So 
we find 



dfAFAA = fB^Apa^d^x = fSjd' 



X 



(B.9) 



where Sj is the rescaled helicity density which we introduce in Sec. 12.41 Finally, we can 
integrate to find the action 

Sj[(l),A^] = I I d^xf^ ^ 



V dtf{^{t))v 



2 J " 2 

And again, the volume average is given by the expectation value 

"^ k^dk 



(B.IO) 



m)) 



V 



{B,A*)6'- 



{2ny 



\A4t,k)\'-\A4t,k)\' . (B.ll) 
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using the results given in Sec. 13.31 

The imphcations for the couphng function from the condition |S'/[0, A^]! < 
|'S'cm[^A']l using the results derived here are discussed in Sec. 14.11 

Appendix C. Asymptotic behaviour of the Coulomb wave functions 

To investigate the late time behaviour of the solution to the mode equation for a power 
law coupling, we need asymptotic expressions for the Coulomb wave function, GqIh, x) 
and FQ{y,x), for x — )■ 0. Based on their expansion in terms of Bessel functions given 
in Abramowitz & Stegun [19], we derive the asymptotes of Gq, Fq, G'q, Fq for x — > 
for arbitrary but fixed y. Actually, we argue that the asymptotic limits for x <^ 2y 
given in 14.6.9-10 (p. 542) of [19] are incorrect. All statements given here were verified 
numerically using Mathematica [59] and Maple [60] . 

We start from the asymptotic expressions in terms of the modified Bessel functions, 
K and J, given in 14.6.7 of [49] for L = 0, and 2y ^ x 



Goiy,x) ^ 2y^Wiy)K^{2y^) (C.l) 

Fo{y,x) ^ -^Wiyr'hi2y^) (C.2) 

where (see Gq given in 14.1.8 of [49] ) 

^ Te-sinhM 1 '^' (C.3) 

ny 

We notice that the approximations given in 14.6.8 of A&S, where 

sinhyry ~ -e'^'^ (C.4) 

is used, are only accurate to better than 1% if y > 1. Since in our case y is arbitrary or 
rather smaller than unity, we cannot employ this approximation. As for the x-derivatives 
of Go and Fq we have verified numerically that the derivatives of the above expressions 
are good approximations, 

G',{y,x) ^ -4yW{y)Ko{2^/2^) (C.5) 

F^(y,x) ^ W{yr'lo{2y^). (C.6) 

Next we use the asymptotic properties of the modified Bessel functions for small 
arguments. As given on p. 375 in [15] 

K,{z) ^ l/z + 0{z), Ko{z) ^ -Hz/2)-^E + 0{z'') (C.7) 

h{z) ^ z/2 + 0{z^), h{z) ~ l + 0{z^) (C.8) 

for z = 2^2xy ^ 1. Using these expressions we find to lowest order in x for x <^ (8?/)"^ 

Go(y,x) ~ W{y), G'^{y,x) ~ 2yW {y)\^{2xy) (C.9) 

Fo(y,x) ^ xlW{y), F'^[y,x) ~ IjW^y). (CIO) 

These asymptotic expressions differ significantly from those given in 14.6.9-10 of [19] 
and it is straightforward to understand why: one arrives at the expressions of 
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when falsely using the asymptotic expansions of the modified Bessel functions for large 
arguments, where Ii{\z\ ^ 1) ~ e^/y/2TTz and -ft'id^z] ^ 1) — e~^^/Tr/{2z). However, 
for X — >■ and fixed y, z = 2y/2xy tends to zero and is not large as assumed for 
approximations 14.6.9-10 of 
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